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We extend the Bell forced dissociation rate model to take account into dynamic disorder. The 
motivation of the present work is from the recent forced dissociation experiments of the adhesive 
receptor-ligand complexes, in which some complexes were found to increase their mean lifetimes 
(catch bonds) when they are stretched by mechanical force, while the force increases beyond some 
thresholds their lifetimes decrease (slip bonds) . Different from our previous model of force modulat- 
ing dynamic disorder, in present work we allow that the projection of force onto the direction from 
the bound to the transition state of complex could be negative. Our quantitative description is based 
on a one-dimension diffusion-assisted reaction model. We find that, although the model can well 
describe the catch-slip transitions observed in the single bond P-selctin glycoprotein ligand 1(PSGL- 
1)— P- and L-selectin forced dissociation experiments, it might be physically unacceptable because 
the model predicts a slip-catch bond transitions when the conformational diffusion coefficient tends 
to zero. 



I. INTRODUCTION 

Traditionally, the forced dissociation rate of noncova- 
lent biological receptor-ligand bonds is described by the 
simple Bell expression m 



where 



fcoff(/) = fcOgexp(/eVfcBT), 



fc°g = fcoexp(-AGVfcBT) 



(1) 



(2) 



is the intrinsic rate constant in the absence of force, /S.G^ 
is the height of the intrinsic energy barrier, ^^ is a projec- 
tion of the distance from the bound state to the energy 
barrier onto external applied force /, fee is the Boltz- 
mann's constant, and T is absolute temperature. Bell 
phenomenologically introduced this expression from the 
kinetic theory of the strength of solids in his seminal pa- 
per more than twenty years ago. Evans and Ritchie J3| 
later put the Bell expression onto a firm theoretical bases: 
the bond rupture was modelled in the framework of 
Kramers theory as thermally assisted escape over one 
or several transition state barriers. The validity of the 
expression has been demonstrated in experiments |3, 14 ■ 
In the Bell forced dissociation rate model, the parame- 
ters including the energy barrier AG''' and the projection 
distance £} are deterministic and time independent. 

On the other hand, we have known that in a 
large variety of chemical and physical areas, associa- 
tion/dissociation processes could be stochastic and time 
dependent. Such processes have been termed as "rate 
processes with dynamical disorder" in an excellent review 
given by Zwanzig |^ in fifteen years ago. The typical 
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examples include cyclization of polymer chains in solu- 
tions [3 , ligands rebinding to heme proteins |7|] , electron 
transfer reactions [g] , electronic relaxation in solutions in 
the absence of an activation barrier |3], and etc. Hence, 
in principle we could expect that dynamic disorder may 
also arise in some forced biological bond rupture exper- 
iments. But it was not until recently that we |10l llll | 
firstly pointed that this disorder might play key role in a 
forced dissociation experiment of the adhesive bond form- 
ing between P-selectin and P-selectin glycoprotein ligand 
1 (PSGL-1) using atomic force microscopy (AFM) [l^ . 
though this experiment was published three years ago. A 
counterintuitive observation has been found: in contrast 
to ordinary biological bonds, the mean lifetimes of the 
PSGL-1— P-selectin bond first increase with initial appli- 
cation of a small force, which were termed "catch" by 
Dembo |l3i | early, and subsequently decrease, which were 
termed "slip". 

The main point of the questions of dynamic disorder 
is that the rate coefficient depends on stochastic con- 
trol variables and is fluctuating in time |3- Correspond- 
ingly, we suggested that the Bell parameters in Eq.^were 
stochastic |ifllll|- As one type of noncovalent bonds, in- 
teractions of adhesive receptor and their ligand is always 
weaker. Moreover, the interface between them has been 
reported to be broad and shallow, such that revealed by 
the crystal structure of the PSGL-1— P-selectin bond [l4l |. 
Therefore it is plausible that the energy barriers AG''- 
or the projection distance £} for the bond are fluctuat- 
ing with time due to either the global conformational 
change or the local conformational change at the inter- 
face. Two possible cases were discussed. In the Gaussian 
stochastic rate model (GSRM) [l3|, we simply assumed 
that ^^ and AG''- had basic properties of random vari- 
ables. We concluded that the catch behavior of some 
biological adhesive bonds could arise from the stronger 
positive correlation between the two stochastic variables. 
Although the model accounted for the catch-slip bond 



transition in a direct and analytic way, it also predicted 
that the mean lifetime of the PSGL-1— P-selectin bond 
was symmetric relative to a critical force fc, while the 
experimental data were clearly skewed towards a large 
force. In addition, the physical justice of the assumption 
of a stationary dissociation process may be doubtful. To 
overcome these shortages, we proposed an alternative and 
more "microscopic" model 11], in which external applied 
force not only presents in the Bell expression, but also 
modulates the distribution of a "inner" conformational 
coordinate. We named it as force modulating dynamic 
disorder (FMDD) model in the following. In contrast 
to the GSRM, we only allowed the energy barrier to be 
stochastic. The quantitative description was based on a 
one-dimension diffusion- reaction equation 7]. We found 
that agreement between our calculation and the data was 
impressive. This model also suggests a new physical ex- 
planation of the catch-slip bond: the transition could 
arise from a competition of the two components of ap- 
plied force along the directions of the dissociation reac- 
tion coordinate and the complex conformational coordi- 
nate; the former accelerates the dissociation by lowering 
the height of the energy barrier, i.e., the slip behavior, 
while the later stabilizes the complex by dragging the 
system to the higher barrier height, i.e., the catch behav- 
ior. 

Even there are apparent advantages of the FMDD 
model, we cannot completely exclude the GSRM, because 
compared to the former, the GSRM does not require that 
force acts on the "inner" conformation coordinate. This 
requirement makes us to conclude that catch-slip bond 
transitions would be altered by the orientation of external 
forces [Ilj. Although such a prediction has a potential 
to account for the bulk experimental observations [I3 
that the optimal binding of P-selectin is critically depen- 
dent on the relative orientations of its ligand, there is no 
such evidence in the exis ting data measured by the sin- 
gle molecular techniques U?, 16]. On the other hand, we 
have viewed the conformational coordinate to be molec- 
ular extension which is coupled to external applied force. 
We also have no strong experimental evidence to sup- 
port this assumption. The inner coordinate has been 
thought to be experimentally unaccessible. In addition, 
the FMDD model did not investigate the possibility of 
a fluctuating projection distance. In the present work, 
we try to give a model that force does not modulate the 
conformational coordinate of the complex while present- 
ing catch-slip behaviors simultaneously. Our model still 
bases on a diffusion-reaction equation. Different from the 
FMDD model, the projective distance £^^ is allowed to 
be fluctuating and negative. Our calculation shows that 
this new Bell rate model with dynamic disorder not only 
derives the same expression of the mean lifetime which 
has been obtained in the GSRM before if both AG-'' and 
^■f are simple linear functions of the conformational co- 
ordinate, but also well fits to the experimental data if 
a slightly complicated piecewise function with two seg- 
ments is used. An unexpected consequence of this model 



is that the transition from catch to slip converts into 
slip-catch transition when the conformational diffusion 
is "frozen" by lowering the diffusion coefficient to zero. 
Interestingly, the FMDD model does not predict such a 
conversion. Therefore, our discussion here would be use- 
ful to design new experiments to distinguish which model 
is physically correct. 

The organization of the paper is as follows. In Sec II, 
we describe the physical picture of the Bell rate model 
with dynamic disorder and give the essential mathematic 
derivations. According to the value of the diffusion co- 
efficient, we distinguish three cases: D -^ cx), D ^ 0, 
and the intermediate D. Analytical solutions can be ob- 
tained for the first two cases, while for the last one we 
resort to numerical approach. In the following section, 
we study two extended Bell models at the three different 
coefficients, in which the barrier height and projection 
distance are linear and piecewise functions of the confor- 
mational coordinate, respectively. A comparison between 
model and the experiments is also performed. Finally we 
give our conclusion in Sec IV. 



II. BELL RATE MODEL WITH DYNAMIC 
DISORDER 

The physical picture of our model for the forced disso- 
ciation of the receptor-ligand bonds is very similar with 
the small ligand binding to heme proteins J]: there is a 
energy surface for the dissociation which dependents on 
both the dissociation reaction coordinate of the receptor- 
ligand bond and the conformational coordinate x of com- 
plex, while the later is perpendicular to the former; for 
each complex conformation x there is a different disso- 
ciation rate constant which obeys the Bell expression 
kosix,f). Higher temperature or larger diffusivity (low 
viscosities) allows x variation within the complex to take 
place. Conversely, at very low temperatures (or high vis- 
cosities) the coordinate is frozen so that the complex is 
ruptured with the rate constant kof{{x, /) without chang- 
ing their x value during dissociation. 



A. Rapid diffusion 

The constant force mode. There are two types of 
force loa ding modes. First we consider the constant force 
case [12, [l3 ■ A diffusion equation in the presence of a 



coordinate dependent reaction is given by 



pr( 

3 



dp{x,t) T^9^P 
dt dx^ 



D d 



k^T dx \ dx 



dV 



P- 



koifix,f)p, (3) 



where p(x, t) is the probability density for finding a value 
x at time t, at initial time p{x, 0) is thought to be thermal 
equilibrium under a potential V{x) without reaction, and 
_D is a constant diffusion coefficient. The motion is under 
influence of a potential V{x) and a coordinate dependent 



Bell expression kos{x, /). Because the above equation in 
fact is almost same with that proposed by Agmon and 
Hopfield except that the rate constant is controlled by 
applied force /, we only present essential mathematic 
derivations. More detailed discussion about the diffusion- 
reaction equation could be found in their original work. 
Following Agmon and Hopfield |3], we substitute 



p{x, t) = Nq(1>{x, t) exp [-V{x)/2kBT] 



(4) 



into Eq.|31 one can convert the Eq.|31into the Schrodinger- 
like presentation, 



dt 



D 



dx"^ 



Uf{x)<t> = -nfi4>), 



(5) 



where A*o is the normalization constant of the density 
function at i = 0, and the "effective" potential 



Uf{x) = U{x) + kos{xJ) 



(6) 



D 



2kBT 



1 



dV 



2kBT \ dx 



dx"^ 



+ kos{x,f). 



We define U{x) for it is independent of the force. In prin- 
ciple Eq. could be solved by eigenvalue technique ■ 
Particularly, at larger D only the smallest eigenvalue 
Ao(/) mainly contributes to the eigenvalue expansion. 
For any given koffix,f), such as the exponential form 
in the Bell expression, there is no analytical Ao(/); in- 
stead a perturbation approach has to be applied .18j . If 
the eigenfunctions and eigenvalues of the "unperturbed" 
Schrodinger operator. 



dx"^ 



^ = -D^ + Uix), 



(7) 



in the absence of ko{f{x, /) have been known to be 

-^4:1 - ->^n4'l (8) 

and the Bell rate is adequately small (compared to the 
diffusion coefficient D), the first eigenfunction (po{f) and 
eigenvalue Ao(/) of the operator Tif then are given by 

Mf) = # + 0[,'H/) + --- (9) 

J<j)'^{x)koff{xJ)(f>l^{x)dxi^ 






and 



>(o) 



^(1)^ 



-^0 ~ ^m 



,(2), 



Ao(/) = X^^" +Xl,'> if) + X'^> {/) + ■■■ 

= A!]+ Ul{x)k,s{xJ)4>l{x)dx + 



(10) 



E 



{Jc^°o{x)k,s{xJ)cl)^{x)dxy 

\0 — \o 
'^0 ^m 



-h 



respectively. Because the system is assumed to be ther- 
mal equilibrium at t = 0, the first eigenvalue Aq must 
vanish. On the other hand, considering that 



0°(a;)(xexp[-T/(a;)/2fcBr], 



(11) 



and the square of 0q is just the equilibrium Boltzmann 
distribution Poq (3^) under the potential V{x), we rewrite 
the first correction of Ao(/) as 



■^0 (./) = Pcqix)koe{x,f)dx, 
Pcqix) ex exp [-V{x)/kBT] . 



(12) 
(13) 



Substituting the above formulaes to Eq.^l the probability 
density function then is approximated to be 



p{x, t) w A^o exp 



2kBT 



exphAo(/)i]<^o(/). (14) 



The quantity measured in the constant force experi- 
ments usually is the mean lifetime of the bond. 



(r)(/) 



t-^dt = 
dt 



Q{t)dt, 



(15) 



where the survival probability Q{t) is related to the prob- 
ability density as 



Q{t) 



p{x,t)dx 



exp 



(a«(/) 



-t- 



(16) 



Therefore the decay of Q{t) is almost a monoexponential 
in large D limits. 

The dynamic force mode. In addition to the con- 
stant force mode, force could be dynamic, e.^., force in- 
creasing with a constant loading rate in the biomembrane 
force probe (BFP) experiment [l^. This scenario should 
be more complicated than that for the constant force 
case. We still use Eq. O to describe the bond dissocia- 
tions by a dynamic force but the constant force therein 
is replaced by a time-dependent function ft. The ini- 
tial condition is the same as before. We firstly convert 
the diffusion-reaction equation into the Schrodinger-like 
presentation again, 



dt 



Hf 



{n + kosix,ft))(f>. 



(17) 



Please note that now we face a time-dependent 
Schodinger operator, Hf^ . 

We know that solving a Schodinger equation with a 
time-dependent Hamiltonian is difficult. The most sim- 
plest situation is to use an adiabatic approximation anal- 
ogous to what is done in quantum mechanics. Hence we 
immediately have [Tg 



(f){x, t) « exp 



{X^{ft.) + B{t'))dt' 



where the "Berry phase" , 

B{t) = I Mft)§-^Mft)dx, 



Mft), (18) 



(19) 



and 4>o{ft) is the first eigenfunction of Hf^. We are not 
very sure whether the assumption is reasonable or not in 
practice. But one of tests is to see the agreement between 
the calculation and data. We apply the perturbation ap- 
proach again to obtain the eigenvalues and eigenfunctions 
of the time-dependent operator. Hence, we have (f>oift) 
and Ao(/t) by only replacing the Bell expression in Eqs.|Hl 
and II 1)1 with fcoff (2;, ft)- The Berry phase then is approx- 
imated to 



B{ft) 



E 



1 



(j)o{x)kos{x, ft)(l>m{x)dx X 



dk, 



4>"o{x)-^CJ^'rni^)dx. 



dt 



(20) 



Finally, the survival probability for the dynamic force is 
given by 



Q{t) « exp 



{\'i\U') + B{Jt, 



dt' 



(21) 



Different from the constant force mode, the data of 
the dynamic force experiments are typically presented in 
terms of the force histogram, which corresponds to the 
probability density of the dissociation forces p{f) 



Pif) 



dQ 
dt 



'dl 
dt 



(22) 



Particularly, when force is a linear function of time / — 
/o + rt, where r is loading rate, and zero or nonzero of 
/o respectively corresponds to the steady- or jump-ramp 
force mode [l^, we have 



P{L 



exp 



1 r 



X\^\f) + B{f) + - 



(23) 



(\'^i\f)+B{n + ---)df 



B. Slow diffusion 

At lower temperatures (or higher viscosities) the con- 
formational coordinate is, or the diffusion coefhcient 
D ^ Q. In this limit, the solution of Eq. Ofor the con- 
stant force mode is simply given by 



p{x, t) = _Poq(a;) exp [~kos{x, f)t] 



(24) 



Compared to the density at larger diffusivity, the proba- 
bility density in this limit closely depends on the initial 
distribution. Consequently, the survival probability Q{t) 
as an integral oip{x, t) is a multiexponential decay which 
is remarkably distinct from the monoexponential decay 
at larger D limitation (Eq. Ilfci|) . The mean lifetime of 
bond is given as 



(t)(/) = / dxpeqix)kj{xj) 



(25) 



For the dynamic force mode, the probability density is 
very simple, 



p{x, t) = _Pcq(a;) exp 



kos{xJt')dt' 



(26) 



and the probability density of the dissociation force for 
the ramping force at the slower diffusion is 



PifJa) = I dxpcciix)-kosix,f) X 



exp 



r 

/ 

/o 



kos{xJ')df 



(27) 



When the diffusion coefficients are between the two 
limiting cases, in general there are not analytical solu- 
tions to Eq. O We have to employ numerical approach 
(SSDP ver. 2.6) ^. 



III. TWO TYPICAL EXAMPLES 

In the present work, we only consider the bound diffu- 
sions in a harmonic potential 



V{x) 






(28) 



where k^ is the spring constant. Ti. then reduces to a 
harmonic oscillator operator with an effective potential 



U{x) 



-'-^ '^X I '^X'*^ 



2k^T \2kBT 



-1 



Its eigenvalues and eigenfunctions are 



A° = nDk^/kBT 



(29) 



(30) 



and 



^(y) = 2-/\-'/\nl)-'/'e-^'f^H,M, (31) 



respectively, where y = (kx/'^kBTY^^x, and Hn{y) is the 
Hermite polynormials 18]. In principle, we can construct 
various functions of the barrier heights and projection 
distances with respect to the conformational coordinate 
X. In the present work, we only study the two examples, 
the linear functions and the piecewise functions with two 
segments. The former should be the simplest functional 
dependence on the coordinate. While the later is complex 
enough to compare with the real data. 



A. Rapid diffusion limit: D ^ 00 



Linear AG-'- and £^^. In this limit, the higher or- 
der corrections Aq , {n > 2) are negligible, and the sur- 
vival probability of bond, Ea. 1161 is monoexponential de- 
cay. Consequently, the mean lifetime of a bond is simply 



(r)(/) = 1/a[)^^ (/). We start from the case that AG* and 
^■t are hnear functions of the conformational coordinate, 



AG*(x) = AGl + kgX, 
^ix) = Cl + hx. 



(32) 



Substituting the above equations and the harmonic po- 
tential into Eq. El we have 



^(/) = fcoexp 



X exp 






-/3AG +^-^-^ 







fcg A;^ / kx £,Q 



a 



f^c/'^x 



(33) 



where /? = l/k^T. We rewrite the above equation into a 
more compact form, 



Ao(/) = fcoexp 



X exp 



-/3AG, 



(3' 



K„ 



{^ - PKg^f 



f3'K^ ( (JKg^-^o 
2 y PKi 



2Kx 

2' 



(34) 



by defining new variables as following, 



K.^ 



/3fc. 



K — 2_ 

' ~ (3kx 






^9? 



fikx 



(35) 



Eq. 1^ has the same expression with the average dis- 
sociation rate formally obtained in the GSRM (Eq. 6 
in Ref. [13 and x\ instead of ^q used therein). Conse- 
quently, Kg , K^ and Kg^ could be viewed as the variances 
and covariance of and between the two stochastic vari- 
ables AG-I- and (^ at the same time point. Interestingly, 
it also means that we cannot get the absolute value of the 
parameters fc^,, kg and fcj by observing the force depen- 
dence of the dissociation rates. Although the above equa- 
tion is not new, there are still several points deserving to 
be presented. Firstly, the derivation of Eq. |^ does not 
require the forced bond dissociation process stationary as 
the GSRM assumed. Then the correlation coefhcient of 
the two stochastic variables is only 1 or —1 according to 
whether kg and fc^ have the same signs or not, whereas 
in the GSRM the coefficient is arbitrary. It is expected 
because both AG-'- and ^ here are the linear functions of 
the same coordinate. Finally, ^q could be minus. We arc 
not ready to discuss the general behaviors of Eq. |^ in 
detail because it has been performed in the GSRM. We 
only emphasize that the rate indeed exhibits the catch- 
slip bond transitions, when the following conditions hold: 
both the distance and the barrier height are stochastic 
variables, and ^ = f3Kg^ — ^q > 0. Because the force is 
always positive, the rate then will firstly decreases with 
the increasing of the force, and then increases when the 
force is beyond a threshold, /c = £^l/(3K^. We define 



AG*,-/^* 



AG*o-/^o 




FIG. 1: Schematic diagram of tlie piece-wise functions AGj = 
AG^{x) — f^^{x) witli respect to the coordinate a; at a given 
force /. The symmetric curve represents the thermal equilib- 
rium distribution of the conformational coordinate under the 
harmonic potential V(x). 



the new parameters £,1 and /c since they have the same 
dimensions of distance and force. 

Following the GSRM, we combine the six parameters 
in Eq. I^into only three: the threshold /c, the variance 
K^, and a prefactor N which is defined as 



N — ko exp 



-l3AGt, + ^Kg-^K^{f-fS^ 



4 



(36) 



Then the mean lifetime has a Gaussian-like function with 
respect to force. 



(r)(/) = iV-^exp 



0^ 



(37) 



We respectively fit the data of the forced dissociations 
of the single PSGL-l-P- and L-selectins [13 bonds. 
The fitting results and the parameters will be showed in 
the next section. Because the Gaussian-like expression is 
symmetric relative to the threshold /c, and it is clearly 
inconsistent with the experimental observations, in the 
following section, we study a slightly complicated case. 

Piece-wise AG-^- and ^-^ functions. Compared to the 
former, the current case should be more useful because it 
describes the experiments better. We assume that both 
the barrier and distance are piecewise function with two 
segments: for x < xt, they follow Eq. [221 otherwise. 



AG*(a;) = AG^ = AG*(a;&), 



— ft — ft/ 



e{x) = e,=eixb) 



(38) 



Fig.^is the schematic diagram of the function AG-'-(a;) — 
f£,^{x) at a given force. Substituting the functions into 



Eq. El we have 
ko 



Ao(/) 



■ exp 



X exp 



-(3AGl 



Pkl 




2h_ 
kf/kx 



pkx ( k^kg I 
' so 



2fc| 



erfc 
ko 



fJ'^X I ^^ kg if 

Z \ kx rZ'T. 



■ exp 



/3AG +/3/e 



X erfc 




where the complementary error function is 



erfc(x) 



e ''^ dx. 



(40) 



Ahhough the above rate expression seems to be more 
complex than that in the case of the linear functions, it is 
really a simple weighted sum of Eq. 1331 and the standard 
Bell expression, Eq. ^ This function would yield the 
major characteristics of the catch-slip bond transitions 
if we assumed that fc^ > 0, ^^ > 0, and kxxl :^ 0: (i) 
if force is smaller. Eg . IS^ reduces to Eg. IS51 because the 
asymptotic behaviors of the error function are 



erfc(— oo) = 2, erfc(+oo) = 0, 



(41) 



respectively; (ii) if force is larger enough, according the 
asymptotic expansion of the error function for large pos- 
itive X, 



erfc(a;) 



1 

2^ 



(42) 



it is easy to prove that Eq. 1351 tends to the standard 
Bell expression with a positive projection distance. This 
behavior has been observed in experiment |2J early. To 
fit the data, we rewrite Eq. [351 into 



Ao(/) 



Ko X { 

exp 



li'K, 



^f! 



X exp 



xerfc 



exp 



2Xb 



f3^Kg 



(3'K^ 



0^ 
2 

/ 



(/ - fc? 



P^Kg 



P^K,J, + 2X,'I^^' 



by introducing additional definitions, 




Y /^^^ ^0 ( n\r<t 

Xb = Xb\l—, Ko = yexp [-fj AG}, 



(43) 



where we have replaced AGl and ^ by the linear func- 
tions. We see that we cannot determine parameters kx, 
kg and fc^ by fitting the forced dissociation experiments. 
Totally there are five parameters presenting in the cur- 
rent model, iiT^, Kg, Xb, Ko, and /c- In practice, we fur- 
ther reduce the model to the simplest case by assuming 
Xb = 0. The fitting curves for the forced dissociations of 
PSGL-1— P-selectin and -L-selectin complex are showed 
in Fig. [21 We see that the agreement between the model 
and the data is quite good. 
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FIG. 2: The mean lifetime vs. force for the bonds of dimeric 
P-selectin and L-selectin with the monomeric sPSGL-1 (the 
square) lia| and the dimeric PSGL-1 (the circle). The reader 
is reminded that the data for the dimeric PSGL-1— P-selectin 
bonds has been rescaled to single bond case |20j. The solid 
lines are the results of Eg. 1391 where the parameters are as fol- 
lows: for P(L)-selectin case, kq « 0.09 (4.60) s"\ Kg ^ 18.5 
(3.39) fcB^^ Ki- « 2.0 (0.06) nm^ and /c = 11.0 (27.7) pN. 
While the dashed lines are from the Gaussian-like Eq. 1341 
where the parameters are as follows: for P(L)-selectin case, 
iV ^ 0.76 (0.085) s-\ K(_ « 1.0 (0.03) nm^ and /c -^ 12.3 
(40.0) pN. Compared to the mean lifetime vs. force in the 
FMDD model, the present curve decays slowly after the tran- 
sition. 



More challenging experiments to our theory are the 
force steady- and jump-ramp modes J16l|. In 13 — > oo 
limit, Eq.j^ni reduces to 



ui f f \ Ao(/) 
Pyj, Jo) = exp 



1 



/ 



/o 



Mf')df' 



(44) 



According Eq. ^^ the probability density of dissociation 
force is uniquely determined by the rate Ao(/). Because 
both the present work and FMDD model describe the 
same group of experiments, here we are not ready to 
discuss the general characteristics of P(/, /o); they have 
been done in FMDD Tl"!. When we test Eg. 1^ with 
the parameters obtained from the constant force data, it 
apparently deviates from the dynamic force data [ig • It 
must emphasize that the experiments of the constant and 
dynamic force modes were performed by the same exper- 
imental group and the same single bond sPSGL-1— P- 
selectin. Similar problem has been met in our previous 



FMDD model and two-pathway and one well model [23| ■ 
We [ij argued and suggested that the experiment of 
dynamic force modes really measured the histogram of 
the dissociation force of the dimeric bond PSGL-1— P- 
selectin complex instead of single bond case: the two in- 
dependent bonds share the same force and fail randomly. 
We still make use of this assumption in the present work. 
It is easy to relate the probability density of the dissoci- 
ation force for the dimeric bonds, Pd{f, /o), to the single 
one by 



PdifJo) = Pif/2Jo/2r 



(45) 



Fig. is the finally result. We see the theoretical pre- 
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FIG. 3: The probability density of the dissociation forces 
Pd{f,fo) under the different loading rates predicted by our 
theory (solid curves) for the PSGL-1— P-selectin complex 
(symbols) 16]. The reader is reminded that each curve re- 
quires an adjustable normalization parameter to fit each his- 
togram. 



diction is satisfactory. Of course, further experiments 
are needed to prove the correctness of this dimeric bonds 
assumption. 



B. Intermediate diffusivity: numerical approach 

Even if the Bell rate model with dynamic disorder can 
well describe the experimental data in quantitative and 
qualitative ways, it does not mean that our model is real; 
Our theory might only work in a specific experimental 
condition, e.g., in the larger D limit. We should test the 
model in a broad range of experimental conditions. Al- 
though until now there are not new experiments about 
the catch-slip bond transitions, we theoretically study the 
general behaviors of Eq. |31 In this and next sections, we 
consider the Bell rate model in the presence of dynamic 
disorder at intermediate and zero D, respectively. In 



principle, the existing experiments would be easily mod- 
ified to test our prediction, such as by decreasing exper- 
imental temperature, increasing viscosity of solvent, and 
etc. Different from the limiting cases, we solve Eq. O 
for intermediate diffusion coefficients D via a numerical 
computing software (SSDP ver.2.6) [13 • Here we focus 
on the piecewise barrier and distance case because. We 
give the reasonable parameters, k =0.09 s~^, k^ =1.0 pN 
nm-\ kg =8.61 pN, k^ =0.70, ^l = 0.64 nm, and Xb = 
nm, which can recover the parameters applied to fit the 
data of PSGL-1— P-selectin bond. The initial distribu- 
tion is set equal to the equilibrium distribution with the 
harmonic potential Eq. |2H1 

We calculate p{x, t) for two diffusion coefficients D and 
two forces; see Fig. 0] When the diffusion coefficient is 
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FIG. 4: Probability density p{x, t) at five times for two forces 
and diffusion coefficients: D — 10^^" and 10^^^ cm^/s for 
the first and second row, and / = 10 and 25 pN for the first 
and second column, respectively. The times are: 0, 175, 937, 
1090, and 2700 ms. The distribution at 2700 ms in the panel 
B are too small to be showed. 



small, the diffusion along x is practically frozen. Hence 
p(a;, t) will decay more on the side that k(x) is larger. 
Because both /S.G^ and £} are piecewise functions with 
two segments, we easily estimate the force, /e = kg/k^, at 
which all decay rates are independent of the coordinate x. 
Because f=10 pN and 25 pN are respectively smaller and 
larger than /o = 12.3 pN given the current parameters, 
for the smaller force, the p(x, t) decays more on the left 
side, while for the larger force, p{x, t) decays more on the 
right side. Fig.0p and D indeed indicate this prediction. 
In contrast, if the diffusion coefficient is large, the dissoci- 
ation is slow compared to process of the conformational 
diffusion. The thermal equilibrium distribution of the 
conformations should be maintained during the courses 
of the dissociation "reaction" . It is cleanly seen in Fig.^JV 
and B. Because the single-molecule forced dissociation 
experiments usually measure the bond survival probabil- 



ity, we calculate Q{t) for the same two diffusion coeffi- 
cients and five forces; see Fig. |31 For the larger diffusion 



and for the piecewise barrier and projection distance case, 
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FIG. 5: Survival probabilities vs. time at five forces for the 
two diffusion coefficients. 



coefficient (Fig|Sp), Q{t)s are almost monoexponential 
decays with time. We see that the decay first slows down 
with force increasing initially, and then speeds up when 
the force is beyond a threshold. These properties have 
been discussed in the rapid diffusion limit. More inter- 
esting case is from the smaller diffusion coefficients, e.g., 
D = 10~^^ cm^s~^ in Fig. |SJI), where the main charac- 
teristics are completely different from the larger D case. 
First, the survival probabilities exhibit highly stretched 
multiexponential decays. The trend becomes more ap- 
parent at larger forces. Second, the decay first speeds 
up and then slows down with force increasing. This ab- 
normal behavior surprises us. It possibly means that the 
catch-slip bond transition in rapid diffusion limit converts 
into slip-catch transition when the diffusion is very slow 
along the conformational coordinate. In order to under- 
stand this behavior better, we consider the nondiffusion 
limit in which Eq. |31has analytic solution. 



(r)(/) = -|-exp 



P^^o + 2 IT - -^ 1 -7— + ^0 



X exp 
erfc 

/Cn 







2 k^^-tAj 



■ exp 



2 



erfc 




(47) 



respectively. We see that, when force is larger enough, 
A:^ > required in I? — *■ oo case leads the mean lifetime 
of the piecewise case in the nondiffusion limit to increase 
rapidly in an exponential way of /^. It proves our esti- 
mation. In addition, for the case of the linear barrier and 
projection distance, this catch behavior at larger forces 
even does not depend on any parameters. Fig. shows 
how the mean lifetime of the PSGL-1— P-selectin bond 
changes with the decrease of D, where same parameters 
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FIG. 6: Mean lifetimes of tfie PSGL-1-P-seIectin bond for 
six diffusion coefficients. The solid and dashed curves are 
respectively calculated by Eqs. 1391 and 1471 in D — > cxd and 
D ->0 hmits. 



C. nondiffusion limit: D — * 

In this limit, according to Eq. [SsJ we get the mean 
lifetimes for the linear barrier and projection distance 
case fEa. 15^ . 



are used. Compared to the counterintuitive catch-slip 
bond transitions, the slip-catch transitions presented here 
is physically inconceivable. 



IV. CONCLUSIONS 



(r)(/) = fc-^exp 



X exp 




Stimulated by the discovery of the catch-slip bond 
transitions in some biological adhesive bonds, several 
theoretical models have been developed to account for 
the intriguing phenomenon. According to the difference 
of the physical pictures, they can be divided into two 
classes. The first class is based on discrete chemical 



reaction schemes. They all assumed that there were two 
pathways and two or one energy well [T^ l20l l2ll | . The 
applied force alters the fractions of the two pathways. 
The other class is based continuum diffusion-reaction 
models and mainly proposed by us (Tfl, UM ■ We cannot 
definitely distinguish which theory or model is the 
most reasonable and more close real situations, because 
their theoretical calculations all agree with the existing 
experimental data. In our opinion, the continuum 
models may be more attractive. In addition that the 
chemical kinetic models are clearly oversimplified, the 
most important cause is that the continuum models 
include an additional parameter, the diffusion coefhcient 
D [ij. The presence of the coefficient obliges the 
continuum diffusion-reaction models to face a strict 
test: it must be physically reasonable in the whole 
range of D. One of examples is that the present Bell 



rate model with dynamic disorder may be physically 
unacceptable because it predicts a slip-catch transition 
when the coefficient tends to zero. We have pointed our 
that the discrete chemical schemes may a mathematical 
approximation of the continuum model [ll|. It is easy 
to see that the two pathway and one energy well model 
with a minus projection distance 2(J| just corresponds 
the present model in the rapid diffusion limit. Therefore 
it is very possible that this chemical kinetic model also 
has a similar flaw. In contrast, FMDD model passes 
this test ^^. Of course, future experiments will finally 
decide which model is correct. 
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